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Abstract 

Cylindrical algebraic decomposition (CAD) is an important tool for the investigation 
of semi-algebraic sets, with applications in algebraic geometry and beyond. We have 
previously reported on an implementation of CAD in Maple which offers the original 
projection and lifting algorithm of Collins along with subsequent improvements. 

Here we report on new functionality: specifically the ability to build cylindrical al¬ 
gebraic sub-decompositions (sub-CADs) where only certain cells are returned. We have 
implemented algorithms to return cells of a prescribed dimensions or higher (layered sub- 
CADs), and an algorithm to return only those cells on which given polynomials are zero 
(variety sub-CADs). These offer substantial savings in output size and computation time. 

The code described and an introductory Maple worksheet / pdf demonstrating the 
full functionality of the package should accompany this report. 

This work is supported by EPSRC grant EP/J003247/1. 

1 Introduction 

This report concerns ProjectionCAD: a Maple package for cylindrical algebraic decomposi¬ 
tion (CAD) developed at the University of Bath. The extended abstract [18] at ICMS 2014 
describes how this package utilises recent CAD work in the RegularChains Library of Maple, 
while still following the classical projection and lifting framework for CAD construction. The 
present report is to accompany the release of ProjectionCAD version 3, describing the new 
functionality this introduced. The report should be accompanied by the code described and 
an introductory Maple worksheet / pdf demonstrating the full functionality of the package. 
The previous two versions of ProjectionCAD are hosted alongside similar reports documenting 
their functionality [niiiTi. 

Version 3 introduces functionality for cylindrical algebraic sub-decompositions (sub- 
CADs): subsets of CADs sufficient to describe the solutions of a given formulae. Two distinct 
types are provided, whose theory was developed in m- The first type contains only those 
cells of a certain dimension and higher, reducing both the output size and computational 
time by giving only output of the required generality. We have implemented both a direct 
and recursive algorithm to build these layered sub-CADs. The second type contains only those 
cells on which given equations are satisfied (lie on a prescribed variety). When building a 
CAD for a formula with an equational constraint then only these cells can contain the solution 
set. These variety sub-CADs clearly reduce the output, and can also reduce computation time 
depending on the rank of the variety relative to the variable ordering. 
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We continue the introduction by reminding the reader of the necessary background mate¬ 
rial on CAD and summarising our previous work with ProjectionCAD. The following sections 
then describe the new functionality in ProjectionCAD version 3. 

1.1 Background on CAD 

A cylindrical algebraic decomposition (CAD) is: a decomposition of meaning a 
collection of non-intersecting cells whose union is M”; algebraic meaning each cell can be 
described by a finite sequence of polynomial relations; cylindrical meaning that with respect 
to a given variable ordering, the projection of any two cells onto a lower dimensional space 
(with respect to the ordering) is either equal or disjoint. CAD was introduced by Collins [11] 
who provided an algorithm which given a set of polynomials would produce a CAD which was 
sign-invariant: so each polynomial had constant sign in each cell. Collins developed CAD 
as a tool for quantifier elimination in real closed fields but CAD has since found applications 
as diverse as robot motion planning [281 etc] and algebraic simplification [131 etc.]. 

Collins’ algorithm (see for example [I]) has two phases. The first, projection, applies a 
projection operator repeatedly to a set of polynomials, each time producing another set 
in one fewer variables. Together these contain the projection polynomials. The second 
phase, lifting, builds the CAD incrementally. First M is decomposed into cells which are 
points and intervals corresponding to the real roots of the univariate polynomials. Then 
is decomposed by repeating the process over each cell using the bivariate polynomials at a 
sample point of the cell. The output for each cell consists of sections of polynomials (where 
a polynomial vanishes) and sectors (the regions between these). Together these form the 
stack over the cell, and taking the union of these stacks gives the CAD of M^. This process 
is repeated until a CAD of M” is produced. Collins’ original projection operator was defined 
so that the CAD of M” produced using sample points in this way could be concluded sign- 
invariance. The key tool to conclude this was polynomials being delineable in a cell, meaning 
the zero set of individual polynomial are disjoint sections and the zero set of the sections from 
different polynomials are identical or disjoint. More efficient projection operators have since 
been developed to conclude both sign-invariance and other invariance conditions. 

The output of a CAD algorithm depends on the ordering of the variables. In this paper we 
usually work with ordered variables ® = xi (so we first project with respect to Xn 

and continue until we have univariate polynomials in xi). The main variable of a polynomial, 
mvar(/), is the greatest variable present with respect to the ordering. 

All cells are equipped with a sample point and a cell index. The index is an n-tuple of 
positive integers that corresponds to the location of the cell relative to the rest of the CAD. 
Cells are numbered in each stack during the lifting stage (from negative to positive), with 
sectors having odd numbers and sections having even numbers. Therefore the dimension of a 
given cell can be easily determined from its index: simply the number of odd indices in the 
n-tuple. Note that the sub-CADs we discuss later are index-consistent, meaning a cell in a 
sub-CAD will have the same index as if it was produced in the associated complete CAD. 

Since Collins published the original algorithm there has been much research into improve¬ 
ments with a summary of the developments over the first twenty years given by [12]. Key 
developments since then include the use of certified numerics [261 [19] in the lifting phase, pro¬ 
jection operators for studying multiple formulae [3] and building CADs by first decomposing 
complex space and refining to real space, instead of projecting and lifting [iniiaiiKis]. 
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1.2 The ProjectionCAD package 

ProjectionCAD is a Maple package developed at The University of Bath to implement CAD 
via projection and lifting. Its name was chosen to distinguish it from the CAD algorithm 
in m which is distributed with Maple as part of the RegularChains library [24112UI etc.]. 
Nevertheless, the package makes much use of procedures from the RegularChains Library 
and the motivation and details for this are given in [18| . 

In [16] we described our implementation of both McCallum’s and Collin’s algorithms to 
produce sign-invariant CADs. In particular, we highlighted that this was the only imple¬ 
mentation to offer order-invariant CADs (meaning each polynomial has constant order of 
vanishing on each cell) and delineating polynomials [22[ [7| which modify the lifting phase 
to allow McCallum’s projection operator to be applied more widely. The latter meant that 
some examples of unnecessary failure in Qepcad [6] could be avoided, while the former was 
necessary for the extensions in the second release. These were described in m and allowed 
for CADs invariant with respect to an equational constraint or CADs which were truth-table 
invariant (TTICADs) with respect to a list of formulae (meaning each formulae has constant 
Boolean truth value on each cell). The TTICAD theory was developed in (S] S] and is based 
on an reduced projection operator analogous to McCallum’s reduced operator for equational 
constraints [23]. In addition it was noted that the reduced projection theory allows also for 
improvement to the lifting phase, leading to ProjectionCAD offering lower cell counts than 
Qepcad even for examples with only one equational constraint. The second release of the 
code also included heuristics to help with choices of problem formulation following [III [5]. 

2 Layered sub-CADs 

Define cells in a CAD with the same dimension as a layer and let I be an integer, 1 < ^ < n-|-l. 
Then an ^-layered sub-CAD is the subset of a CAD of M"' with cells of dimension n — i 
for <i < i. They were introduced by the present authors in |27j (although there had been 
previous work where only cells of full-dimension were computed mm)- 

Algorithm [1] describes how an Ulayered sub-CAD may be produced with the main idea 
to run the projection phase as normal (step[T|), but then truncate in the lifting phase when 
cells of a high enough dimension cannot be produced. The algorithm can run with different 
projection operators: that of Collins [T] or McCallum [22] to build sign-invariant CADs for 
polynomials (as detailed in m); or operators to build truth-invariant CADs for formulae 
using equational constraints [23| and truth-table invariance nn (as detailed in ini)- 

We build the CAD of the real line as normal (step [2|) but then for each successive lift we 
first check the dimension of the cell to be lifted over (step[6|) and only proceed to generate the 
stack if cells of dimension I or higher in M"’ could be produced from it. The stack generation 
command is detailed in |18j and makes use of the RegularChains stack generation code 
requiring a preprocessing step. Some projection operators can incur theoretical failure for 
input that is not well-oriented, in which case this is identified during stack generation and 
FAIL returned. A final loop is used (steps [14] — [T7|l to remove any lower dimensional cells 
that were produced as part of a stack. The approach was verified in m- 

As previously described in [16] there are various output formats available. At the least 
each cell is represented by an index (positioning it in the CAD) and a sample point (encoded 
as a regular chain and bounding box) while an algebraic description is also available. Of 
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Algorithm 1: LCAD(F, a;,^,ProjDp): Algorithm to produce Alayered sub-CADs. 
Input : A choice of projection operator Pro j Dp; input F{x) (polynomials or 
formulae) in ordered variables X] and an integer ^ E (l,n -|- 1). 
Output: An Alayered sub-CAD for F, or FAIL if F not well-oriented. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


Run the projection phase and set P[i] to be the projection polynomials with mvar x*; 
Set P[l] to be the CAD of obtained by isolating the roots of P[l]; 
for i = 2,..., n do 
Initialise T>[i\ := [ ]; 
for c E T>[i — 1] do 

Evaluate dim := EaGc.index (« “od 2); 
if dim > i — I — 1 then 

S := GenerateStack(P[i], c); 
if S =FAIL then 

I return FAIL; // Input is not well oriented 

else 

1^ Add the cells in S to V[i]] 


13 

14 

15 

16 
17 


Initialise P := [ ]; 
for c E P[n] do 

Evaluate dim := Easc.index (« 
if dim > n — i then 
Add c to V] 


mod 2); 


18 return P; 


particular note it the intuitive piecewise construction (described in [8]) which highlights the 
tree-like nature of CAD. This has been adapted in ProjectionCAD to display layered sub- 
CADs with the truncated branches clearly indicated. 

A simple example of using LOAD in ProjectionCAD now follows, appearing as it would 
in a Maple worksheet, except that the sample points have been replaced by SP for brevity. 


> LCAD([x''2+y''2-l] , 1, [y,x] 
( SP 


method=McCallum, output=piecewise); 

X < —1 

branch = truncated 




SP 

y < —V— + 1 


* * * 

branch = truncated 


SP 

—V—x^ + 1 < y < V—x^ -P 1 

— 1 < X < 1 

* * * 

branch = truncated 


SP 

\/-x2 + 1 <y 


* * 


branch = truncated 


SP l<x 


Here a sign-invariant sub-CAD for the defining polynomial of the unit circle has been 
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produced. The number 1 as the second input indicates that only one layer of cells is to be 
produced; those five cells with full dimension. Two simplifications have been in this sub- 
CAD compared to a complete CAD. First when lifting over the CAD of the real line the two 
points {x = ±1) have not been lifted over: they failed the dimension check at step [6] and were 
then excluded from the output during steps [H] — [T71 Second, those cells in the stack over 
the interval (—1,1) on the real line withut full dimensional were also discarded in the final 
loop. Note that while the second simplihcation reduced the output, the hrst also reduced the 
computation time since no real root isolation was required over 2 of the 5 cells in M^. 

We could have instead asked for two layers of cells to receive the output below. 


]. 

2, [y,x], method=McCallum, output=piecewise); 

'SP 


X < —1 


'SP 

y < 0 


< 

* * * 

branch = truncated 

X = —1 


SP 

V 

0 <y 



'SP 

y < —\/—x^ P 1 



SP 

y = -^-x^ P 1 


< 

SP 

— \/—x‘^ pi <y < y/—x^ P 1 

— 1 < X < 1 


SP 

y = pV-x^ P 1 



SP 

\ 

\/—x^ -t-1 < y 



'SP 

y < 0 


< 

* * * 

branch = truncated 

X = 1 


SP 

V. 

0 < y 


SP 

V 


1 < X 

missing 

cells are the isolated points x 

= ±l,y = 0. These would be 


included in a 3-layered sub-CAD, which is itself a complete CAD. Hence in this case LOAD 
would give exactly the same output as CADFull (described in [T6]l. 

Note that each of LOAD command is distinct: re-running the projection phase and recom¬ 
puting the higher dimensional cells. ProjectionCAD also contains a recursive command to 
build layers incrementally. It takes additional input as a list of cells already computed and 
sections that were previously discarded (in the final loop). It then lifts over those discarded 
sections to produce cells of one dimension higher, as well as further discarded sections. This 
recursive algorithm is described fully and verified in m- 

Algorithm [T] allows for projection operators other than McCallum’s. In particular, it can 
be used with the TTICAD projection operator of [3] to build truth-table invariant sub-CADs, 
as for the following problem. 


> fl := x*2+y*2+z*2-l: gl:=x*y*z-l/4: 

> f2 := (x-l)*2+(y-l)*2+(z-l)*2-l: g2:=(x-l)*(y-l)*(z-l)-l/4: 

> PHI := [ [fl, [g2]], [f2, [g2]] ]: 

The syntax here means that PHI dehnes two formulae: the hrst with EC dehned by /i and 
additional constraints dehned also by gi ; and the second with EC dehned by /2 and additional 
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constraints defined by also g 2 - The TTICAD command could produce a decomposition sign- 
invariant for fi and /2 and also for each gi where the corresponding /* is zero. 

> TTICAD(PHI, [z,y,x]); 


497 

The LTTICAD command can be used to build only the prescribed layers of cells: 

> LTTICADCPHI, 1, [z,y,x]), LTTICAD(PHI, 2, [z,y,x]), 

LTTICADCPHI, 3, [z,y,x]), LTTICAD(PHI, 4, [z,y,x]); 

93, 299, 455, 497 

So for example, producing only the full-dimensional cells gives less than a fifth of the output. 
It also takes less than half the computation time. 

3 Variety sub-CADs 

An equational constraint (EC) is an equation, / = 0, logically implied by the truth of a 
Tarski formula. They may be explicit, like / = 0 A (/>, or implicit as / 1/2 = 0 is in 

(/i = 0 A (pi) V (/2 = 0 A 4>2). 

The presence of an EC can reduce the number of projection polynomials required [23l [3l etc.] 
and ProjectionCAD already has the related commands ECCAD and TTICAD (described in |17j). 
Given a Tarski formula with EC / = 0 we define a variety sub-CAD as a truth-invariant 
sub-CAD for the formula consisting only of cells lying in the variety defined by / = 0. 

Algorithm [2] describes how a variety sub-CAD may be produced. We assume that all 
factors of the EC have the same main variable as the overall system (xn) and so the EC is 
used only for the first projection and final lift. Building variety sub-CADs outside of this 
restriction is considered in [27] but not yet implemented in ProjectionCAD. 

The projection proceeds as normal in stepdj as does lifting to a CAD of in step[2| 
For the final lift we generate stacks only with respect to the EC in step [T] noting that the 
projection theory used should ensure this is sufficient to conclude truth invariance (see |1] for 
details). Further, only the sections in those stacks are retained for the output in step [TTl as 
only these describe the variety. Since the underlying CAD algorithms that utilise ECs can 
return theoretical failure for input that is not well oriented we assume this is tested for during 
stack generation. This approach was verified in m- 

Consider the formula (/? := (x^ -|- — 1 = 0) A {xy — | > 0) A (x^ — > 0). To identify the 

regions in where this is true we could build a sign-invariant CAD for the three polynomials. 

CADFulK [x~2+y~2-l ,x*y-l/4,x~3-y"2] , [x,y] , method=]yicCaIIum) : nops(7o); 

161 

However, since we see the formula has an EC we can use this to build a truth-invariant CAD 
for (p with the ECCAD command: 
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Algorithm 2: VCAD(<^,/, x): Algorithm to produce variety sub-CADs. 

Input : A choice of CAD projection Pro j Dp; input F{x) (polynomials or formulae) in 
ordered variables x] an EC / = 0 (in which all factors have same mvar as F). 
Output: A (truth-invariant) variety sub-CAD for or FAIL if F not well-oriented. 


1 Set P to be the output from applying Pro j Op to E; 

2 Compute a CAD of for P and assign to V] 

3 if V' = FAIL then 

4 1^ return FAIL; // P is not well oriented 


5 

6 

7 

8 
9 


Initialise T> := []; 

for c € P' do 

Calculate S := GenerateStack(/, c); 
if 5 = FAIL then 
1^ return FAIL; 


// E is not well oriented 


10 

11 


if I S'! > 1 then 

1^ Add those cells in S with last entry of their index even to E; 


12 return E; 


ECCAD( [x''2+y''2-l, [x*y-l/4,x~3-y~2]] , [x,y]): nops("/o); 

73 

This more than halves the number of cells. The VCAD command goes further and returns only 
those cells where the EC is satisfied (a sub-CAD) more than halving the output again: 

VCAD( [x~2+y~2-l, [x*y-l/4,x~3-y~2]] , [x,y]): nopsC"/); 

28 

Note that in all 3 cases it will be necessary to then test a sample point from each cell to see 
where the formula is true or false. Although the cells outputted by VCAD guarantee the EC 
is satisfied, the two inequalities are only guaranteed truth-invariant, and so may still be false. 

4 Additional functionality in ProjectionCAD version 3 

Combining sub-CADs 

We can build sub-CADs which are both variety and layered using the LVCAD command. In fact, 
we can combine all the ProjectionCAD theory to build layered variety truth-table invariant 
sub-CADs with the LVTTICAD command. More examples and details are given in m, while 
in |28] such a combination was used to solve a long-standing robot motion planning problem. 

Cell distribution in CADs 

Tools are provided to compare the growth in CADs by layer, that is, the distributions of cells 
in CADs by dimensions. Analysis in [29] suggests a common distribution with little variation 
for the underlying problem, meaning layered sub-CADs can be used as heuristics to guide the 
construction of complete CADs. 
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Using sub-CADs to avoid theoretical failure 

Algorithms based on McCallum’s theory can only run on input that is well-oriented (a 
definition particular to each individual operator, see for example mm)- These conditions 
are checked during lifting and if not satisfied ProjectionCAD gives a warning. For example: 

> f := a*e+b*d+c*e+d+e: vars:=[a,b,c,d,e]: 

> CADFulK [f] , vars, method=McCalluin ): 

Warning, The input is not well-oriented (nullification on cell [1, 2, 2]). 

The output cannot be guaranteed correct. 

We can build a 1- or 2-layered sub-CAD for the polynomial without triggering the warning: 

> LCAD( [f],1, [a,b,c,d,e], method=McCallum,failure=err): 

> LCAD( [f] ,2, [a,b,c,d,e] , method=McCallum,failure=err) : nopsCyo'/o), nops(y); 

48,148 

Hence we can get descriptions of not only the full dimensional cells but also the next dimension 
down without running into well-orientedness issues. 


5 Summary 

We have described the new functionality present in the third release of ProjectionCAD, 
which focusses on algorithms to build sub-CADs. The underlying theory is discussed further 
in [281 [291 [27] and the present report should be accompanied by the code and an introductory 
Maple worksheet demonstrating the full functionality of the package. A key topic for future 
work is building variety sub-CADs when the EC does not have all factors in the main variable. 
Various approaches are laid out in Section 2.1 of [27] and an analysis of these is ongoing. 
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